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Abstract 

We have performed an ecliptic survey of the Kuiper belt, with an areal coverage of 8.9 square degrees to a 



50% limiting magnitude of 



Sloan 



= 24.7, and have detected 88 Kuiper belt objects, roughly half of which 



received follow-up one to two months after detection. Using this survey data alone, we have measured 
the luminosity function of the Kuiper belt, thus avoiding any biases that might come from the inclusion of 
other observations. We have found that the Cold population defined as having inclinations less than 5° has a 
luminosity function slope acoid = 0.82±0.23, and is different from the Hot population, which has inclinations 
greater than 5° and a luminosity function slope aHot = 0.35±0.21. As well, we have found that those objects 
closer than 38 AU have virtually the same luminosity function slope as the Hot population. This result, 
along with similar findings of past surveys demonstrates that the dynamically cold Kuiper belt objects likely 
have a steep size distribution, and are unique from all of the excited populations which have much shallower 
distributions. This suggests that the dynamically excited population underwent a different accretion history 
and achieved a more evolved state of accretion than the cold population. As well, we discuss the similarities 
of the Cold and Hot populations with the size distributions of other planetesimal populations. We find 
that while the Jupiter family comets and the scattered disk exhibit similar size distributions, a power-law 
extrapolation to small sizes for the scattered disk cannot account for the observed influx of comets. As well, 
we have found that the Jupiter Trojan and Hot populations cannot have originated from the same parent 
popuation, a result that is difficult to reconcile with scattering models similar to the NICE model. We 
conclude that the similarity between the size distributions of the Cold population and the Jupiter Trojan 
population is a striking coincidence. 
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1. Introduction 

The size distribution is one of the most funda- 
mental properties of a planetesimal population. As 
the size of an object is primarily determined from 
its accretion and collisional disruption histories, the 
size distribution can reveal important information 
on the accretion and collisional histor y of that pop- 



ulatio n (for a recent example, see iBottke et al 
2005h . 



Unlike the closer populations such as the aster- 
oid belt, whose proximity has allowed the accu- 
rate measurement of their size distributions (see for 



example Uedicke and Metcalfe! . Il998t Uewitt et al 



"Accepted by Icarus August 2nd, 2010. 
Email address: fraserwOgps.caltech.edu (Wesley C. 
Fraser ) 



2000) , the distance to the Kuiper belt has prevented 
an equally detailed determination of its size distri- 
bution. Observations thus far have demonstrated 
that for the Kuiper belt as a whole, the size distri- 
bution for objects with diameters, D > 200 km, 
is well described by a power-law. At some size 
50 < D < 150 km, the size distribution rol l s-over 
to a shallower distribution (jBernstein et all 120041 
Fuentes et ail 120091: iFraser and Kavelaarsl l2009h 
Smaller than the roll-over, the size distribution ap- 
pear s to remain shallow to ob jects as small as D ~ 1 
km (jSchlichting et al.l . l2009h . 



While the accuracy of current measurements pre- 
vents a detailed modelling of the history of ob- 
jects in the region, some insight has already been 
gained. The steepness of the large object size 
distribution implies that for this population, ac- 
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cretion was a short-lived p henomenon, like l y not 



more than a f e w 100 Myr (Gladm an et all 12001 



Kenvonl l2002t IFraser et all l2008h . The paucity 



of the observed belt, and the existence of the 
largest known members demonstrates that the belt 
has undergone significant mass depletion, losing 
as mu ch or more than 99% of its primordial 
mass (IStern and Colweli 19971 IKenvon and Luu . 



1998[ IJewitt and SheppardL l2002t iFuentes et al 
2009). The existence of the roll-over at sizes 
larger than D > 50 km suggests that the 
belt has undergone significant collisional com- 
minution in a region of significantly increased 



density compared to today (IKenvon and Bromle 



2004 
2009). 



Bcna videz and Campo Bagatin, 2009; Frascr 



As well, recent observations have suggested that 
the size distribution of those dynamically cold 
(low inclinations and eccentricities) Kuiper belt ob- 
jects is different than that of the dynamically hot 
(large inclinations and e c centricities) population 
|Levison and Sternl . l200lh . IBernstein et all (|2004l ) 
found that the size distribution of objects with in- 
clinations, i < 5° was steepe r than the size distri- 
bution of objects with i > 5°. Fuentes and Holmanl 
(2008) found similar results showing that the cold 
population had a steeper slope than the mixed pop- 
ulations as a whole. This result implies that the hot 
and cold populations are genetically separate pop- 
ulations, that have significantly different accretion 
and evolution h istories. 

The results of IBernstein et al 1 (<2004h were drawn 
from samples of objects compiled from many dif- 
ferent surveys. This practice was necessitated by 
available data. All surveys had either shallow lim- 
iting magnitudes and a large areal coverage, or 
vice-versa. The result was that most surveys did 
not have a sufficient range of objects for which 
the size distribution could be accurately measured 
from that survey alone. The practice of using 
data from multiple surveys opens the results to the 
possibility of being affected by calibration issues 
and variations in sky density which could lead to 
an incorrect measurement of the size distribution. 
When these effects were properly account e d for , 
Fraser et all (|2008l ): IFraser and Kavelaarsl (|2009l ) 
found that these data could not be used to reliably 
test differences in subpopulations of the Kuiper 
belt. 

Here we present the results of a new Kuiper belt 
survey. By virtue of the survey's design, a large 
number of objects were discovered and followed to 



determine their inclinations, over a range of sizes 
sufficient to measure the size distribution for the 
hot and cold population without the use of other 
surveys. In Section [2] we present our observations, 
and our data reductions and discovery techniques. 
In Section[3]we present the analysis of our results, in 
Section 0] we discuss the implications of our results, 
and we end with concluding remarks in Section [5l 



2. Observations and Reductions 

2.1. Discovery Observations 

The discovery observations were made on Octo- 
ber 18, 2009 (UT) with Suprime-cam on the 8.2 m 
Subaru telescope (jMivazaki et all 120021 ). Suprime- 
cam is 10-chip mosaic camera, with a field-of-view 
of roughly 34'x27' with ~ 15' chip-gaps, and has 
a pixel scale of 0.2". The observations consisted 
imaging 34 fields, each visited three times in the 
r' filter with the camera long axis oriented horizon- 
tally in RA, with each visit consisting of a single 200 
s exposure. The first and last images of a field were 
separated by roughly 6.5 hours. The observed fields 
were all within ~ 1.5° of the ecliptic and within 
~ 10° of opposition at the discovery epoch. The 
total areal coverage of our survey after accounting 
for chip-gaps was 8.93 square degrees. Details of 
the individual fields are shown in Table [T] and Fig- 
ure[TJ The observations included multiple images of 
the Dl and D4 fields of the Canada- France-Hawaii 



Teles cope Supernova Legacy Survey (|Astier et al 



2006) at various airmasses to provide both photo- 
metric and astrometric calibrations of the discovery 
images. 

The images were pre-processed with standard 
techniques; the bias levels were removed, and a mas- 
ter bias frame, set as the average of 10 bias images 
was removed from all science frames. A master sky- 
flat was produced from the science frames using a 
clipped median filter, and was removed from the 
science images producing background variations no 
larger than ~ 1% across a mosaic. 

Photometric calibrations were done on a chip-by- 
chip basis using the Mega-pipe so urce catalog of 
the Dl and D4 fields (|Gwvnl . [2008h . The zeropoint 
was found to be identical for all chips within un- 
certainties, with a value, Z rgubaru — 27.50 ± 0.02. 
The conversion from the Suprime-cam instrumental 
magnitudes to those of Mega-pipe was found to be 
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m{r')suprime = m(r') M ega - 0.018(g' Mega - r' Mega ). 

(1) 

Atmospheric extinction was measured from the sci- 
ence frames. The airmass correction was found to 
be 0.08 magnitudes per unit increase in airmass. 
Comparison of common background sources in each 
triplet revealed no large seeing or transparency vari- 
ations throughout the night (see Figures [2] and [3]) ; 
the night was photometric. 

The Suprime-cam field of view is highly distorted, 
especially near field edges, requiring calibration to 
ensure accurate astrometr y of moving sources in 
the observations. Scamp (Bertin, 2006) was used 
to measure second order spatial distortions of the 
Suprime-cam field, from the Dl and D4 fields on a 
chip-by-chip basis. Little variance in the distortions 
was found over the airmass range of these observa- 
tions. Thus, a master distortion map for each chip 
was produced from all images of the Dl and D4 
fields. Absolute calibration of each science image 
was done using the USNO-B catalog, and the mas- 
ter distortion map, resulting in residuals of roughly 
0.3" . Relative astrometry between the three images 
of each field - a triplet - was performed by matching 
sources between an image, and the reference image, 
chosen as the image of the triplet with the low- 
est absolute residuals. This resulted in an excellent 
common astromctric solution for the triplet, with 
residuals of < 0.06" between all three images. This 
calibrations procedure ensured reliable astrometry 
to both maximize the discovery efficiency of moving 
sources, as well as secure the possibility of follow-up 
on future dates. 

2.2. Moving Object Search 

The Suprime-cam data were utilized to search for 
moving objects in the Solar system, in and beyond 
the orbits of the gas-giant planets. To characterize 
the detection efficiency of our moving object detec- 
tion method, artificial moving point sources were 
implanted in the observations. 

The stellar point-spread function (PSF) was gen- 
erated on an image-by-image and chip-by-chip basis 
from 10-15 hand-selected, visually inspected point 
sources for each chip. It was found that a spatially 
constant PSF was sufficient to model the stellar 
shape across a chip. The PSFs were created using 
the tools in the doaphot package of IRAF and for 
each image consisted of an average Moffat profile 
with look-up table of the PSF stars of that chip. 



Random artificial sources were generated on Sun- 
bound orbits, with semi-major axes, 18 < a < 1000 
AU, eccentricities, e < 0.6, and inclinations, i < 
90°. Anomaly, periapse, and nodal angles were 
chosen at random to ensure that the sources fell 
in the images and could be either approaching or 
receding from their nodes. 20-30 artificial sources 
were planted in each chip with apparent magnitudes 
in the range 21 < r 1 < 26.5 with image-to-image 
flux variations matching those measured from 30-40 
bright stars common to each image of a triplet, thus 
ensuring that any extinction variations throughout 
the observations were accounted for. If an artifi- 
cial source's motion would cause it to drift more 
than 0.2 pixels during an exposure, it was broken 
up into multiple fainter sources with 0.2 pixel spac- 
ing and a total flux equal to that of the artificial 
source. Thus, trailing effects were fully accounted 
for in our search. 

Moving objects were identified by their motions 
between ima g es of a triplet. Using Sextractor 
( Bertin et all |2002[ ). all sources in an image were 
tabulated. Stationary sources common to each im- 
age in a triplet were flagged, and ignored from fur- 
ther consideration, leaving only those sources who 
were either moving, or had highly variable fluxes. 
The remaining sources of an image triplet were then 
searched for motion consistent with Sun-bound ob- 
jects. Candidate sources were chosen as any three 
individual point-source detections, one per image, 
whose centroid moved at least one pixel (0.2") and 
at most ~5" between each image. These candi- 
date sources were then further fil tered by fit_radec 
( Bernstein and Khushalanl 2000l) . If the best-fit 
orbit to the candidate had a chi-squared larger than 
7.5, the candidate was rejected. The results of this 
were approximately 16000 candidate sources which 
were then visual inspected by an operator for final 
acceptance or rejection. 

The majority of visually inspected candidates 
were false, primarily caused by cosmic ray impacts, 
or extended sources which were identified by Sex- 
tractor as more than one source in some images. 
Of the 16,000 candidates, roughly a third were 
matched with implanted artificial sources. From 
these sources, the performance of the search rou- 
tine was characterized. The results are shown in 
Figure HI As can be seen, the detection efficiency is 
well described by the familiar functional form 



1 — tanh 



(2) 
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where A, m*, and g are the maximum efficiency, 
half-maximum magnitude, and half-width param- 
eters. The survey achieved a 50% limiting magni- 
tude for discovery of r'=24.7. The search was found 
to be sensitive to objects as distant as ~ 900 AU, 
beyond which, insufficient motion was exhibited for 
reliable detection. No significant variation in search 
efficiency with inclination or eccentricity was found 
(see Figures [5] and [6J . 

The search resulted in the detection of 88 sources. 
The sources' distances and inclinations were de- 
termined with fit_radec. Photometry was mea- 
sured using standard aperture photometry tech- 
niques with corrections for airmass applied (see Fig- 
ure [3]) . Aperture corrections were measured from 
the bright stars from which the PSFs were gener- 
ated. Photometric uncertainties were determined 
from the function 



Ar' = 6 + 7IO 1 



(3) 



where 5 is a constant representing the extinction 
stability of the night in question, 7 is a constant 
proportional to the parameters of the detector (see 



Fraser et all |2008, for a derivation), and Z is the 



telescope photometric zeropoint. This function was 
fit to the difference in measured and planted mag- 
nitudes of the artificial sources. The best-fit pa- 
rameters are (£,7) = (0.028,3.25). The best-fit S is 
in agreement with the scatter in extinction exhib- 
ited in Figure El Along with the measured fluxes 
with uncertainties sampled per measurement from 
Equation [3j these results are presented in Table [2 
One detected object was identified as previously 
known Plutino 2004 VT75 which was not known 
to be in the observed fields prior to detection. This 
detection confirmed the accuracy of both the photo- 
metric and astrometric calibrations. The detected 
position of 2004 VT75 was ~ 0.2" from that pre- 
dicted from the Minor Planet Center. As well, this 
source was detected across different chips demon- 
strating that the astrometric fits for different chips 
of the same mosaic image are reliable in an absolute 
sense. 

2.3. Follow-up 

For a subset of the objects, follow-up obser- 
vations were performed with the Low Resolution 
Imaging Spectrometer (LRIS) on the Keck-I tele- 
scope on the nights of November 21st, 22nd and 
December 21st, 2009 (UT). LRIS is a dual channel 
(blue and red) camera, with each camera, a 2 CCD 



mosaic having approximately a 6X8' field of view, 
with a ~ 25" gap, and a pixel scale of 0.135". This 
particular instrument is convenient for our follow- 
up observations, as the positional uncertainty el- 
lipse of our targets at the time of follow-up was 
slightly smaller than the areal coverage of a single 
detector CCD. 

The follow-up observations consisted of pairs of 
images at the positions of each target, predicted 
from fit_radec, with temporal spacing between each 
image sufficient to detect the object's motion. By 
virtue of the dual channel imaging abilities of the 
camera, images in both the g and r LRIS filters were 
taken simultaneously. Exposure times were scaled 
from the flux of each source in the discovery images, 
and were chosen such that each source would have 
a photometric signal-to-noise ratio of ~ 10 in each 
image. Images of the Dl and D4 Supernova fields 
were taken for calibration purposes. 

Due to time constraints, only those objects who's 
distance might place them beyond 39 AU, and 
with magnitude in discovery images brightward of 
r' = 24.5 were followed. 41 objects satisfied these 
constraints and were pursued in the follow-up ob- 
servations. 

Due to the geometry of the follow-up observa- 
tions, asteroid confusion was a concern. The pre- 
dicted rate of motion of the follow-up targets was 
used as a diagnostic of whether or not a detected 
moving object was the correct source. The second 
CCD of the LRIS detector provided a measure of 
this asteroid confusion rate, as no Kuiper belt ob- 
jects were expected to fall within this region of the 
images. In the November follow-up data, 3 aster- 
oids with brightnesses and rates of motion consis- 
tent with our Kuiper belt targets were found in the 
second CCD for all follow-up exposures. With the 
exception of two follow-up pointings, our objects 
were not confused, as only one moving target per 
field was identified with a rate of motion consistent 
with the targeted source. This is consistent with 
our asteroid confusion rate determined from the 
other CCD. Special attention for the two objects 
with confused follow-up was made during the De- 
cember observations. From the december data, the 
interloper in the November observations was easily 
rejected achieving reliable links between the discov- 
ery and follow-up images for all follow-up targets. 
No confusion occurred during the December follow- 
up. 

Our follow-up efforts were rewarded with 100% 
success rate for all 41 targets satisfying our dis- 
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tance and magnitude cuts. The conditions on the 
November nights were photometric and of moder- 
ate seeing, with full- width at half maximum < 1.5 . 
On the night in December, while it was photomet- 
ric, the seeing conditions were mediocre, reducing 
the observing efficiency compared to the November 
nights. As such, all of the 41 objects for which we 
attempted follow-up have month-long arcs, while 
only a small subset of 10 have arcs longer than this. 
The results of the follow-up efforts are presented in 
Tabled 

Astrometric and photometric calibrations of the 
LRIS images were performed using the same tech- 
niques as applied to the Suprime-cam images. The 
astrometric distortions of the LRIS field were small, 
but noticeable, and varied with airmass. Distortion 
maps extracted from the calibration images at air- 
masses similar to the science images provided resid- 
ual errors of ~ 0.3-0.4" with respect to the USNO- 
B catalog. 

For the LRIS blue-channel, the photometric zero- 



point was found to be Z„ 



27.91 ±0.03 with an 



airmass term A g — 0.15±0.05. During the observa- 
tions, the LRIS red-channel detectors were plagued 
by electronics difficulties causing charge transfer 
inefficiencies that were dependent on source flux. 
This effect was found to be as much as ~ 10% in 
engineering images. As such, we adopt a 10% un- 
certainty on the zeropoint of this detector, and find, 
Z rKcak = 28.55 ±0.1. The red-airmass term was 
found to be A r = 0.08 ± 0.03. 

The conversions between the LRIS r and g mag- 
nitudes and the r' and g' Mega-pipe magnitudes 
were found to be 



m(r) Keck = m(r' M ) - 0.2(g' M - r' M ) (4) 



and 



m{g) K eck = m{g' Mega ) + 1.2(g' Mega - r' Mega ). (5) 

The relatively large correction between gKeck and 
g' Mega is caused by the LRIS beam-splitter configu- 
ration which, in the configuration we utilized, trun- 
cated the red side of the gKeck filter while preserv- 
ing the throughput across the rKeck filter. 

2.4- Color Conversion 

To facilitate the combination of the discovery 
and follow-up images, all photometry needed to be 



converted to a sing le system. The Sloan system 
(jSmith et all 120021 ) was chosen. The conversion 
from Megaprime magnitudes to the Sloan system 
are 



^Mega ^ ' Sloan 



0.024Q4 o 



' ' Sloan 



) (6) 



and 



9Mega ~ 9 Sloan C '.153(fl ' Sloan r 'sioan)- CO 

From Equations!]] 21 and[S]and the conversions pre- 
sented above, we find that the conversions from the 
Subaru and Keck photometry to the Sloan system 
are 

r 'sioan = r Keck - O.Q94(g Keck - r Keck ), (8) 



9sic 



gKeck - 0.445(g_R- 

eck TKe 



(9) 



and 



r Slo 



rSubaru ~ 0-019(g K eck ~ r K eck)- (10) 



Using Equations ISHTUI all photometry have been 
converted to the Sloan system. The photometry 
presented in Table [5] are the weighted averages of all 
measurements presented here. When no {gKeck — 
rKeck) colour was available, the average color of all 
those objects that did receive follow-up, (gKeck — 
rKeck) = 1-21 was used to convert to the Sloan 
survey. As the range of Keck colors for our sample 
is 0.68 < (gKeck — rKeck) < 1.96, we estimate the 
error in the conversion caused by using this average 
color is ~ 0.02 magnitudes, negligibly small. 

In Figure E| we present the (g' sloan - r' sloan ) color 
versus inclination for those objects with follow-up 
observations. From this figure, it can be seen that 
our sample exhibits similar behaviour to the Kuiper 
be lt objects in th e Mino r Planet Center analyzed 
bv lPeixinho et al~ ( 20081) . That is, there is a corre- 
lation with colour and inclination; low inclination 
objects are typically redder than higher inclination 
objects. As both our colours and inclinations are 
significantly mo re inaccurate than thos e in the sam- 
ple analyzed by iPeixinho et all (|2008l ). a repeat of 
the same analysis on our data is not warranted, 
as no more insight would be gained than by visu- 
al ly comparing our sam ple to theirs (see Figure 1 
of IPeixinho et all (l2008l)1. 
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3. Analysis 



As found by iBernstein et ah (2004), their seems 
to be some variation in the Kuiper belt size distribu- 
tion with orbital excitation. More specifically, they 
defined two populations, the Cold sample, which 
contains those objects with heliocentric distances, 
38 < d < 55 AU and inclinations, i < 5°, and the 
Excited population as all those Kuiper belt objects 
in the same distance range, but with large inclina- 
tions. IBernstein et al. ( 2004f) found that for bright 
objects (R < 24), the Cold population had a much 
steeper size distribution than the Excited popula- 
tion. This result was found from a sample of Kuiper 
belt objects compiled from a large number of dif- 
ferent surveys, providing a result possibly affected 
b y calibration errors and sky density variations. 

lElliot et al.l (|2005l) analyzed the results of the 
Deep Ecliptic Survey, and found similar results. 
That is, that the cold classical objects have a 
steeper size distribution than the resonant, or ex- 
cited populations. This result however, was drawn 
from a survey without a calibrated detection effi- 
ci ency, making the resu l t untr ustworthy. 

iFuentes and Holman d2008f) per formed a similar 
analysis as Bernstein et al.l ( 2004 ) . but with addi- 
tional survey data. Their findings were similar, ie., 
that the hot population size distribution is shal- 
low er than that of t h e cold population. But much 
like Bernstein et al. (l2004h . they considered data 
from multiple surveys and may have been affected 
by the same calibration errors and sky density vari- 
ations. 



Fraser et al 



Fraser and Kavelaars 



reanalyz e d the survey data considered in 
Bernstein et al. (l2004h along with new observa- 
tions, and introduced a technique which accounted 
for possible errors induced from combining the 
results of different surveys when measuring the size 
distribution. They found that the apparent differ- 
ences between the hot and cold populations could 
not be disentangled from possible effects caused 
by the combination of different surveys. When 
correctly accounting for these effects, no statistical 
difference between the hot and cold population 
was found. Population differences as large as those 
inferred by IBernstein et ah ( 2004 ) might exist, but 
their presence could not be discerned from the 
available data sets. 

The survey we present here, has sufficient depth, 
and number of detections to probe the differences 
in size distributions of the hot and cold populations 



alone, without the need for other data. By doing 
so, we can ensure that the biases that may be intro- 
duced from the use of multiple datasets are avoided 
here providing a clean and reliable test for any dif- 
ferences between the hot and cold populations of 
the Kuiper belt. 

3.1. The Samples 



We want to test the findings of IBernstein et al 
(|2004l ). To that end, we define subsets of our de- 



tections, which are si milar to those considered by 
Bernstein et al. I (l2004h . They are defined as follows: 



1. We define TNO sample, as all the detections 
in our sample. 

2. Like Bernstein et al. (2004) we define the Cold 
population as those objects with 38 < d < 55 
AU and i < 5°. We also further refine this, and 
define the Coldp sample as only those sources 
from the Cold sample who received follow-up 
observations in addition to the discovery ob- 
servations. 

3. We define the HotF and Hot samples as those 
objects with 38 < d < 55 AU and i > 5°, who 
have and have not received follow-up observa- 
tions respectively. 

4. We define the Close population, as those ob- 
jects with 30 < d < 38. We do not make any 
distinction about follow-up, as the majority of 
the Close population was not observed beyond 
the discovery observations. 

We note here that the exact inclination division 
between the Cold and Hot populations is uncer- 
tain. Measurements of the inclination distributions 
of the Hot and Cold samples suggest that the vast 
majority of cold classic al objects wil l have inclina- 



5° (Brown, 200 f; Gulbis et al 



tions less than 

120101 ) . Analysis of the optical colour distribution of 
Kuiper belt objec ts reveals a cor r elatio n of colour 
with inclination (Peix inho et al. . 2008. for a re- 



cent example, see) that suggests there is a divi- 
sion between low-inclination red objects, and low- 
inclination neutral objects at i ~ 12°. Given the 
uncertainty of the inclination division that sepa- 
rates the Cold and Hot populations, we consider 
the full range of inclination divisions testable by 
our observations and primarily discuss the histori- 
cal division of i = 5° unless otherwise stated. 

The Close population is a mix of objects in mean- 
motion resonances with Neptune, the Centaurs, and 
the scattered disk objects. On the other hand, 
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the Hot population consists of hot classical, reso- 
nant, and scattered disk objects. While the exact 
fractional mix of each of these populations cannot 
be determined from our observations, the Minor 
Planet Center can guide us as to which populations 
dominate these subsamples. Analysis of the objects 
in the Minor Planet Center, near opposition in late 
October, within a few degrees of the ecliptic, reveals 
13 objects satisfying 30 < d < 38 AU, 6 of which are 
Plutinos, and 36 objects satisfying d > 38 AU and 
i > 5°, of which only a few are Plutinos, with the 
majority being hot classical and scattered disk ob- 
jects. As the relative mixes of each dynamical class 
are different for the Hot and Close subsamples, any 
differences detected in the size distributions of the 
subsamples can provide insight into the size distri- 
butions of the underlying populations. 

3.2. The Methods 

As the Kuiper belt objects we detect are not re- 
solved, we have no measure of their sizes, and there- 
fore cannot measure the size distribution of our 
sample directly. Rather, we must use an in-direct 
technique, and infer the size distribution from the 
luminosity function. Brightward of R ~ 25, the 
luminosity function is well represented by a power- 
law. That is, the differential surface density of ob- 
jects with magnitude m per square-degree is given 
by 



E(m) = aln(10)10' 



a{m — rn 



.) 



(ii) 



where a is the logarithmic slope of the power- 
law, and m is the magnitude at which there is 
one object per square degree with that magni- 
tude, or brighter. Recent measurements suggest 
that, for the TNO sample, a = 0.73 - 0. 76 and 

" 20091: 



23.4 in the R-band (Fuentesetal 



Fraser and Kavelaarsl 20091 ). It is worth noting 
that the luminosity function has a break, or tran- 
sition from the steep slope for large objects, to a 
shallower distribution for fainter sources. The tran- 
sition occurs around R ~ 25, which is at the tail 
end of our survey sensitivity. As such, the break 
will have little to no effect on our observations, and 
we ignore it in our analysis. 

It can be shown that if the size distribution 
of a population is a power-law, with slope q, ie., 
tjt cx r~ q , then the luminosity function of that pop- 
ulation is given by Eauation llll with slope a = . 
As such, the slopes of the size distributions of each 



population can b e determined from their luminos- 
ity functions (see Fraser et al. , 20081 for a thorough 
discussion of this technique) . 

The fits of Equation [11] is done with a Bayesian 
maximum likelihood technique. Specifically, the 
likelihood 



L{{m}\a,m ) = e~ n * n(™)X(m\ a ,m )d m 
x j^Jcr(mi)E(mi|a,m ) 

(12) 

is maximized over the luminosity function slope and 
normalization parameters, a and m . Here Q is 
the areal coverage of our survey, {to} is the set of 
magnitudes for our discovered objects, with rrn the 
magnitude of object i. o~{mi) is a functional repre- 
sentation of the uncertainty in the magnitude TOj. 
As standard practice, we adopt a gaussian repre- 
sentation, with widths equal to the uncertainty in 
the observed magnitude of each source. 

To avoid potential errors induced at low detec- 
tion efficiencies, we consider only those sources 
with probability of detection greater than 50% and 
truncate our detec tion efficiency below this point. 
Fraser et al.1 ( 20081 ) has demonstrated that truncat- 
ing the efficiency at 50% as we do here causes the 
best-fit slope a to appear roughly 3-5% steeper than 
in reality. This effect however, is much smaller than 
the uncertainty in the best-fit slopes we present. 
For the CoUf and Hotp samples, our follow-up is 
limited to those sources with r' Subaru < 24.5. Thus 
r' = 24.5 is our magnitude cut for the CoUf and 
Hotp samples. In Equation [TT]we utilize only the 
magnitudes from discovery observations alone, as 
utilizing the data from other filters could possibly 
introduce biases into the results. We note however, 
that when utilizing all available flux measurements, 
the change in the results was significantly smaller 
than the uncertainties in the best-fit parameters. 

To test the quality of the best-fits, we utilize the 
Anderson-Darling statistic, 



(S(m) -P(to)) 2 
P(m) (1 - P(m)) 



dP{m) (13) 



where P(m) is the cumulative probability of detect- 
ing an object with magnitude < to, and S(m) is the 
cumulative distribution of detections. We calculate 
the probability, P(A > A & s ) of finding a value, A, 
larger than that of the observations, A f, s , given the 
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best fit-parameters by bootstrapping the statistic, 
ie., randomly drawing a subsample of objects from 
the best-fit power-law, with number equal to that 
detected in a particular sample, fitting Equation [TT] 
to that random sample, and computing A. Values 
of P(A > A t, s ) near indicate that the functional 
form is a poor representation of the data. 

3.3. The Results 

The luminosity functions of the TNO , Cold , 
Hot , and Close samples are presented in Figure |8j 
The results of the fits of Equation [TT1 to the various 
subsamples are shown in Figure |H1 and Table [3] 

As can be seen from Figures [5] and the Cold p 
exhibits a much steeper luminosity function, with 
slope acoid = 0.82 ±0.23 than the Hotp population 
which has slope ctHot = 0.35 ± 0.21. The slopes of 
these two populations differ at more than the 1-er 
level. Both of these samples are well described by 
power-laws, as exhibited by their Anderson-Darling 
statistics (see Table [3]) . Similar slopes within the 
uncertainties are found from the Hot and Cold sam- 
ples. But these values are less trustworthy as the 
inaccurate inclinations of those objects which did 
not receive follow-up cause an uncertain amount of 
mixing between the two subsamples. 

In addition, the Close sample exhibits virtu- 
ally the same luminosity function slope, aciose — 
0.40 ± 0.15 as the Hot F ( and Hot ) sample. As 
the Close and Hot subsamples are made up of dif- 
ferent fractions of the Centaur, resonant, hot clas- 
sical, and scattered populations, these results sug- 
gests that all of these populations exhibit equally 
shallow luminosity functions. If one or more of the 
populations had a luminosity function as steep as 
the Cold population, it would likely be detectable 
as a difference in a between the Hot and Close 
populations, the steeper being the subsample with 
the greatest fraction of the steep dynamical pop- 
ulation. This however, is not apparent. Indeed, 
the slope does not change when considering the 
Hot and Close populations together, resulting in 

UHot+Close =0.40 ±0.12. 

The luminosity function of TNO sample has a 
slope ctTNO — 0.56 ±0.1, which is significantly shal- 
lower than the slope a ~ 0.75 found from other 
recent measurements in the same magnitude rang e 



recent me asurements m the same magnitud e range 
( Fuentes et al. , 2009; iFraser and KavelaarsL [2009). 



This surprising result can be understood relatively 
simply. The TNO sample, is a mixture of the Hot , 
and Cold populations, which our results have shown 



have significantly different slopes. As such, we ex- 
pect the TNO sample to have a slope bounded 
by the slopes of the Hot (ctHot ~ 0.4) and Cold 
(ctcoid ~ 0.8) populations. Our observations were 
made at a ecliptic longitude where plutinos, which 
are primarily members of the hot population, are 
preferentially at perihelion (see above). As such, 
we expect a greater fraction of Plutinos detected 
in our observations, than at other longitudes where 
Plutinos are not at perihelion. As the Hot popu- 
lation has a shallower luminosity function than the 
Cold population, we expect to see a flatter luminos- 
ity function for the TNO sample at this longitude 
than at others. This simple reasoning can account 
for the variation in luminosity function slopes seen 



2008). 



in past surveys (see IFraser et al 

Along with the size distribution differences 
between the Hot and Cold populations , the 



colour distribution ( Teeler and Romanishinl . 
Truiillo and BrownL 120021: IPeixinho et al 



and s ee Figure [7]), inclin ation distribution ( B 



20011 ). binary fraction ([Stephens and Noll 120061) 



2000; 



2008 



rown, 



and lack of lar ge objects on dynamically cold orbits 
( Brown . 20081) suggest that these two populations 
are different. There is no clear evidence however, 
for a clean separation in inclination between the 
two. Rather, it is more likely that a smooth gra- 
dient exists, and as such, the historical division of 
idiv = 5° chosen to separate these two populations 
is relatively arbitrary. As such, we measured the 
luminosity function slopes while varying idiv The 
results are presented in Figure 1101 

As can be seen, the Cold and Hot population 
slopes become significantly different for idiv ^ 5°. 
When considering all objects, not just those with 
follow-up, the slope of the Hot population seems to 
become a constant, with ctHot ~ 0.4 for idiv > 6°. 
This suggests that above 6° the fraction of the pop- 
ulation with a steep size distribution is insignificant, 
and the observed objects come almost entirely from 
the population with a shallow size distribution. As 
the uncertainty on a is still quite large, these ob- 
servations however, are insufficient to confirm such 
a hypothesis. 

To test the significance of the apparent difference 
in luminosity functions of the Hot and Cold popu- 
lations, we utilize two different statistical tests. We 
first utilize the Kuiper-variant of the Kolmogorov- 
Smirnov test as a non-parametric measure of the 
significance i n the observe d difference between two 
populations ( PressL 120021 ). We randomly boot- 
strapped from the cumulative luminosity function 
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of the Cold sample, a sample of objects equal in 
number to the Hot sample. We then calculated 
the KS-statistic of this random sample compared to 
the Cold sample. This process was repeated, and 
from this the probability of finding a simulated KS- 
statistic as large or larger than that found from the 
actual Hot population was determined. We found 
that the two populations were drawn from separate 
parent populations at the 80-90% significance levels 
for inclination divisions idm > 4°. 

While the KS-test provides a non-parametric 
test, a better measure of the significance differ- 
ence of the Hot and Cold populations can be found 
with the knowledge that the two populations are 
well described by power-laws. To that end, we also 
utilized the Anderson- D arling statist i c, as d one in 



the luminosity function of the TNO population ex- 
hibits a roll-over at magn itudes R > 25 i s wel l 



Bernstein et~aH (|2004 : iFraser et all (|2008l) . Fol 



lowing similar procedures as above, by randomly 
sampling from the best-fit power-law of the Cold p 
sample we determined the probability of finding 
a simulated Anderson-Darling statistic as large or 
larger than that found from the actual Hotp pop- 
ulation as compared to the observed CoUf sam- 
ple. We found that the two populations were 
drawn from separate parent populations at greater 
than the 3-er significance for inclination divisions 

idiv ^ 4 . 



4. Discussion 

4-1- Implications for the Kuiper belt 

The observed luminosity functions imply size dis- 
tribution slopes qcoid = 5.1 ± 1.1, quot = 2.8 ± 1.0, 
and qciose — 3.0 ±0.8 for the Gold , Hot , and Close 
populations. The consistency of the Hot and Close 
populations leaves us to consider the two popula- 
tions together, implying a slope for the combined 
population of qnot+ciose = 3.0 ± 0.6. 



Tho ugh the results presented by IBernstein et al 
(|2004h and lFuentes and Holmanl (|2008D were possi- 



bly affected by sky density variations and calibra- 
tion errors, we find similar conclusions, ie., the Cold 
population size distribution is steeper than the Hot 
population size distribution. Given the similarity 
in these repeated measurements of the Kuiper belt 
luminosity function, the results stand confident. 

While over the magnitude range of our obser- 
vations, both the Hot and Cold populations are 
well described by power-law models, other works 
have put into question whether the power-law be- 
haviour extends to all brightnesses. The fact that 



accepted. The results of IBernstein et al.l (j2004) 



however, suggest that the power-law behaviour of 
the Cold population breaks down for the brightest 
objects. Specifically they find that at the bright- 
end, the Cold population luminosity function has 
a ste eper slope than at fain ter magnitudes. Simi- 
larly, MorbideiESaD (|2009l) suggests that the Hot 
and Cold populations are intimate mixtures of two 
primordial populations, one with a steep size dis- 
tribution, and one with a shallow size distribution. 
In such a scenario, both populations should exhibit 
shallow slopes for the largest objects, then turn-up 
to a steeper slope for smaller objects, before finally 
breaking to a shallow slope at some size where col- 
lisional processing has dominated. To that end, we 
test the power-law behaviour of both the Hot and 
Cold populations. 

For the Cold population, a simple extrapola- 
tion of the best-fit luminosity function suggests 
that the brightest object on the sky in this pop- 
ulation should have a magnitude of R ~ 19.5. 
The entire low-latitude Kuiper belt has been sur- 
vey ed for objects to a l i miting magnitude of R < 
21 (jTruiillo and BrownL l2003h . and the brightest 
known Cold object has R ~ 21.3, nearly a full 
magnitud e fainter. This observ ation confirms the 
results of iBernstein et aL (2004). Namely, that the 
power-law behaviour of the Cold population lumi- 
nosity function cannot extend to the brightest tar- 
gets. Rather the power-law must be truncated. 
The exact behaviour cannot be determined from 
the observations we presen t here. The suggestion 
by iMorbidelli et al.l (|2009h that the cold popula- 
tion size distribution will have a shallow slope for 
the largest objects, and a turn-up to the observed 
steep is excluded by this result. 

Similarly, extrapolating the best-fit power-law lu- 
minosity function of the Hot population suggests 
that this population is well described by a power- 
law to its brightest objects. Indeed, when ex- 
cluding the largest objects which are known to 
have different albedos than the objects observed 
here ( Stansberrv et all 120081 ) the absolute magni- 
tude distribution of the biggest obje cts has a slope 



Morbidclli et al 



2008 



compatible with our observations ([Brown 

" 20081) . To test the exactness of 



a power-law requires combination of other survey 
data sensitive to the brightest Hot members. As 
well this analysis should include modelling of the 
latitude distribution of the Hot population. The 
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fact however, remains that the Hot population has 
a flat size distribution over all observable sizes. 

The behaviour of the Hot population luminosity 
function at the faint end is uncertain; the shallow 
slope of a.Hot ~ 0.35 for R < 25 is compatible with 
the lack of faint detections of the Hot population in 
fainter surveys, eliminating the need for a break at 
the faint end. If however, the Hot population slope 
is as steep as the upper limits of our confidence 
interval, a break at magnitudes fainter than R ~ 26 
is still required. The existence of a break in the Hot 
population will only be confirmed by an off-ecliptic 
survey sensitive to objects with R > 26. 

The turn-up, as proposed by Morbidelli et al.l 
( 20091) would produce an abundance of faint objects 
in our observed Hot population. We can test this 
utilizing the Anderson-Darling statistic. We gen- 
erate a fake Hot subsample from a shallow power- 
law with the observed Hot population slope that 
turns up to the best-fit observed slope of the Cold 
population, acoid = 0.8, at some magnitude rriTu- 
We then fit a power-law to that sample, and cal- 
culate the Anderson-Darling statistic. This is re- 
peated, and the probability of finding a statistic 
worse than the observations is calculated versus the 
value of rriTu- The strength of this test is weak- 
ened by the sparse sampling of high-latitude sources 
in our dataset. We still however, eliminate such 
a turn-up to the slope of the Cold population for 
m TU < 24.2 (r') at the 2-a level and < 23.6 at 
the 3- g level. This is incomp atible with the asser- 
tion of Morbi delli et al. I (|2009l) who suggest that the 
Hot population should have the same slope as the 
Cold population over an absolute magnitude range 
of 6.5 < H < 9, corresponding to r' magnitudes, 
22 < t' < 24.5. The significance of the test results 
over the entire magnitude range however, prevents 
this result from being iron clad. Clearly, the possi- 
bility of a turn-up (or down) must be tested from a 
survey which detects a large sample of high-latitude 
objects. 

The difference in size distribution slopes of the 
Hot and Cold populations suggest very different 
histories for these two groups. Interestingly, the 
shallow slope of the Hot distribution is compati- 
ble with a heavily collisionally processed popula- 
tion, either one that has reached collisional equi- 
librium, or in which the largest objects are frag- 
mcnt s of even larger, disru p ted p ri mordial bod- 
ies (lO'Brien and Greenbere . 20031 : iBottke et al 
2010). For collisional evolution to completely re- 



would require extremely high collision rates, seem 
ingly incompatible with plausible pr otoplanetar 
disk densities, and formation scenar i os ( Sternl . Il99 



Kenvon and Bromlev . 2004 : Fraser . 2009h . Rather 



it is likely that this slope is the result of the accre- 
tionary processes that formed the Hot population. 

Unlike the Hot population, the steep slope of the 
Cold population is entirely incompatible with a col- 
lisionally evolved distribution. This slope must be 
the result of accretionary processes as well. 

The shallower slope of the Hot population, and 
the fact that the largest objects of the Hot pop- 
ulation are larger than those of the Gold popula- 
tion, implies that the Hot population achieved a 
more advan ced stag e of a ccretion than the Cold 
population (jKenvonl . l2002h . The incompatibility 
of their slopes implies that these two populations 
underwent different accretion scenarios. It is possi- 
ble, that the Hot population underwent a longer 
duration of accretion than the Cold population. 
To halt accretion for the Cold population would 
require excitation and mass depletion with some 
yet unseen perturber beyond the outer edge of 
the Kuiper belt. The past existence of one or 
more planetary embryos in the Kuiper belt re- 
gion is a possible source of this excitation, and 
has been proposed to e xplain some of the Kuiper 
belt dynamical featu r es rtGladman and Chan . 20061 : 



Lvkawka and Muk ai. 2008). 



shape the size distribution of the Hot population 



Another possibility is that the Hot population 
underwent accretion in a more dense region of the 
protoplanetary nebula. This would result in more 
rapid accretion than compared to the Cold pop- 
ulation, allowing the Hot population to grow to 
larger sizes before the process was halted. Under 
this scenario, the Hot and Cold populations could 
then have similar accretion timescales. As currently 
favoured formation scenarios suggest the Hot pop- 
ulation was scattered from a regi on closer to the 
Sun than their current locations (jMalhotraL [1993; 
iGomed . l2003 iLevison et all l2008h , where the pro- 
toplanetary disk might have been more dense, this 
scenario seems likely. 

However the Hot and Cold populations came 
about, their incompatible size distributions imply 
different accretion scenarios. Once accretion was 
finished, these objects were excited and emplaced 
onto their current orbits, creating the architecture 
of the current Kuiper belt. Given the dynamically 
cold nature of the Cold population, it seems plau- 
sible that this population formed in-situ. What 
ever the mechanism(s) ultimately responsible for 
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the belt, our results show that little mixing between 
the Hot and Cold populations has occurred. 

4-2. Comparison with Other Populations 

There is a striking similarity in the size distribu- 
tions of the Cold population, and the Jupiter Tro- 
jans. Both have similarly steep slopes, (qcold ~ 
5.1 and qjr ~ 5.5), and both exhibit breaks to 
shallowe r slopes at roughly the same object di- 



amet er (|Jewitt et al.l . 120001 : iFraser and Kavelaarsl 
2009). This result argues strongly against the Tro- 



jan formation mechanism of the so-called NICE 
model in which the Trojan populations and the 
hot and cold Kuiper belt populations are all scat- 
tered into their current re gions from the same pri- 



2005 



mordial disk population j Morbidel li et al 
Nesvornv and Vokrouhlickvi 2009h. If this see 



nario were true, iMorbidelli et al. r fcoogh has shown 
that over the magnitude range of our observations, 
the hot population should have the same slope as 
that observed for large Jupiter trojans. Using the 
Anderson-Darling statistic, we tested the probabil- 
ity of the Hot sample being drawn from a luminosity 
function with the slope of the Jupiter trojan lumi- 
nosity function, ajr — 0.9. We found that the Hot 
sample could not be drawn from the Jupiter trojan 
population at greater than the 3-er significance. As 
the size distributions of these two populations in 
the size range considered have not evolved signif- 
icantly s ince they were em p laced in their current 
regio ns ( Davis et al. . 2002 : Fraser and Kavelaarsl 
2009h . we must conclude that the hot and Jupiter 
trojan populations must have different progenator 
populations. This result and the lack of mixing be- 
tween the hot and cold populations is difficult to 
reconcile with the NICE model. It seems likely the 
hot, cold, and trojan populations formed by sepa- 
rate means. If this is true, the similarities between 
the cold and trojan populations are quite a coinci- 
dence. 

Another similarity is seen between the Hot pop- 
ulation (quot ~ 3) and the Jupiter family comets, 
who exhibit a size distribution slope of qjFC ^2.8 
albeit over a smaller size range, 1 S, D S, 10 km 



(jTancredi et all 120061: IWeissman et~"alll2009h . This 



resemblance is interesting as one likely source of 
Jupiter family comets are scattered disk objects - 
members of the Hot population - which have fallen 
into the inner Solar system under gravitational per- 
turbations from the gas-giants, suggesting that the 
scattered disk size distribution might be a power- 



law for D > 1 km. We consider this possibility 
here. 

While the observations we present do not mea- 
sure the size distribution of the scattered disk di- 
rectly, the similarity in slopes between the Hot and 
Close subsamples implies that the scattered disk 
cannot have a size distribution significantly differ- 
ent than that of the Hot sample as a whole. Addi- 
tionally we cannot determine the total number scat- 
tered disk objects from our observations. Rather, 
we turn to other estimates which suggest that, if 
their size distribution is a power-law with slope of 
q ~ 3, then there are roughly 10 7 — 10 8 scattered 



disk objects with 1 < D < 10 km (jTruiillo et al 



2000l : Parker and Kavelaarsl |2010L Schwamb, per- 
sonal communicat ion) . 

Simulations by IVolk and Malhotra ( 2008 ') sug- 
gest that if the scattered disk is the sole source of 
the Jupiter family comets, then there must be at 
least 10 9 objects in that population to account for 
the current flux of Jupiter family comets through 
the inner Solar system, implying that either the 
scattered disk is not the sole source of the Jupiter 
family comets, or that the extrapolation of a power- 
law to f ~ 1 km is unreasonable; it must be that 
the size distribution of scattered disk objects steep- 
ens significantly in the D ~ 10—100 km range. 
Further observations are required before this fea- 
ture can be detected. 



5. Conclusions 

We have performed an ecliptic survey, and have 
detected 88 Kuiper belt objects with magnitudes 
21 < r' 3loan < 25.2. A subset of these objects have 
received additional follow-up observations allowing 
us to accurately determine their inclinations and 
distances. Using these data, we have measured the 
size distribution of the Hot and Cold subsamples, 
historically defined as those objects with inclina- 
tions above and below 5° respectively. This mea- 
surement, which is independent of any previous ob- 
servations has confirmed that the Cold population 
has a much steeper luminosity function, with slope 
Oicoid = 0.82 ±0.23, than the Hot population, with 
slope oiHot = 0.35 ± 0.21. 

The observed luminosity functions imply differ- 
ent size distributions for the Hot and Cold pop- 
ulations. The size distribution slopes of the the 
two subsamples are qcold = 5.1 ± 1.1 and quot = 
2.8 ± 1.0. 
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In addition, we have found that the Close popu- 
lation, which is defined as those objects with helio- 
centric distance, d < 38 have the same luminosity 
function slope as the Hot subsample, demonstrat- 
ing that these two share a similar size distribution. 
In addition, our findings suggest that the dynamical 
populations which make up both the Hot and Close 
populations must all have similar size distributions. 

The primary consequence of these findings is that 
the Cold population, which consists primarily of 
cold classical Kuiper belt objects, has a separate, 
and distinct accretion history from the Hot popula- 
tion, requiring either potentially different accretion 
timescales for the two populations, or formation 
in different locations of the protoplanetary disk. 
These observations reveal the similarities in the size 
distributions between the subsamples of the Kuiper 
belt and the Jupiter family comets and Trojans, 
suggesting a connection between the formation and 
subsequent evolution of the small body populations 
of the outer Solar system. 
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Table 1: Suprime-cam Discovery field centers. Coordinates are 
presented in the J2000 epoch. 
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Woijoo SDSSs 


R (AU) 2 


i (°) 2 
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25.2 ± 0.1 


25.0 ± 0.1 




42.0 ±3.0 


10.0 ± 9.0 
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24.55 ± 0.07 


24.52 ± 0.07 




42.0 ±3.0 


2.0 ± 8.0 
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24.75 ± 0.08 


24.72 ± 0.08 




40.0 ±3.0 


20.0 ± 10.0 


Nobj_3 
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obj_79 


24.85 ±0.09 


24.78 ±0.09 




50.0 ±3.0 


1.0 ±6.0 


obj_80 


24.58 ±0.07 


24.54 ±0.07 




42.0 ±3.0 


5.0 ±8.0 


obj_81 


24.03 ±0.05 


24.01 ±0.05 




33.0 ±2.0 


6.0 ±6.0 


Nobj_82 


24.19 ±0.05 


24.13 ±0.05 


0.78 ±0.09 


43.1 ±0.1 


0.31 ±0.01 


obj_83 


23.06 ±0.03 


22.95 ±0.03 




35.0 ±3.0 


30.0 ± 10.0 


obj_84 


24.98 ±0.1 


24.94 ±0.1 




41.0 ±3.0 


2.0 ±8.0 


obj_85 


24.54 ±0.07 


24.52 ±0.07 




30.0 ±3.0 


17.0 ±9.0 


Nobj_86 


23.87 ±0.04 


23.84 ±0.04 


0.68 ±0.09 


42.0 ±3.0 


2.1 ±0.6 


obj_87 


25.1 ±0.1 


25.1 ±0.1 




45.0 ±3.0 


8.0 ±9.0 



Tabic 2: Discovery Subaru magnitudes and Sloan magnitudes. 



1 Objects prefaced with N and D have follow-up observations in Novembe r and December respective ly. 
2 Barycentric Distance and ecliptic inclination determined from fit_radec (Bernstein and Khushalanill200(j|) . 
3 Object distance and inclination taken from the Minor Planet Center. 
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Table 3: Best-fit power-law parameters. 



Sample 


a 


TOo(r') 


P(A > A obs ) 


TNO 


0.56 ±0.1 


22.9 ±0.4 


0.64 


Coldp 


0.82 ±0.23 


23.8 ±0.3 


0.28 


Cold 


0.79 ±0.20 


23.8 ±0.3 


0.47 


Hotp 


0.35 ±0.21 


24.3 ±0.7 


0.75 


Hot 


0.35 ±0.19 


24.3 ±0.7 


0.65 


Close 


0.40 ±0.15 


23.6 ±0.6 


0.68 



Table 4: The slopes of the TNO , Hot , Cold , and Close samples (rows 3, 5, and 6) are likely a few percent steeper than in 
actuality. This bias however, is much smaller than the uncertainties of those parameters (see Section 13.211 . 



Fig. Captions 

Fig. [JJ Diagram of the fields observed for moving objects during the Subaru observations. These fields 
were chosen to avoid bright stars. The field center coordinates are presented in Table [TJ 

Fig [2] Point source full- width at half maximum versus hour angle of all Subaru discovery observations. 

FigEl Extinction versus airmass for the Subaru discovery observations. Each point is the average difference 
in magnitude of bright stars between common fields; each of the 35 discovery triplets has two data points. 
The line is the best-fit linear extinction law, with slope 0.08 magnitudes per unit airmass. These data 
demonstrate the photometric nature of the discovery observations. 

Fig |U Detection efficiency versus artificial source r' Subaru magnitude. Error-bars are 1-sigma poisson 
limits. Best-fit parameters of the efficiency function given by Equation [5] are presented at the top of the 
figure. 

Fig [5l Detection efficiency versus artificial object inclination. From top to bottom, all planted sources, 
and those with r'< 23, 23 < r' < 24, 24 < r' < 26. This figure demonstrates that, given a particular object 
magnitude, our survey was equally sensitive to all inclinations. 

Fig |51 Detection efficiency versus artificial object heliocentric distance. From top to bottom, all planted 
sources, and those with r'< 23, 23 < r' < 24, 24 < r' < 26. This figure demonstrates that, given a particular 
object magnitude, no significant variation of detection efficiency with distance is apparent. Exceptions are 
at the extrema of the survey; our search was only sensitive to objects with apparent motions 0.2 — 10" /hr 
roughly corresponding to distances 15 < r < 800 AU. 

Fig [7] Sloan colour versus inclination for those objects which received follow-up observations. 

Fig El The differential luminosity functions of the various samples defined in Section [3. II offset in units of 
2 for clarity. Diamonds: TNO . Squares: Coldp . Triangles: Hotp . Circles: Close . Error-bars are the 1-cr 
poisson limits for the number of objects in each bin. Lines are the best-fit power-laws, with values given in 
Tabled 

Fig [3 The 1 (solid), 2 (dashed) and 3-ct (dotted) likelihood contours of the fits presented in Figure E] and 
Table [3] The corresponding subsample is labeled in the top right of each panel. 

Fig [10] Best- fit slopes a of the Cold (squares) and Hot (circles) populations versus the inclination division 
idiv separating the two. Left is using only those objects with follow-up. Right is using all observations. The 
1-er range in otciose is bracketed by the two horizontal lines. 
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